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ABSTRACT 

For a molecular cloud clump to form stars some transport of magnetic flux is required from the 
denser, inner regions to the outer regions of the cloud, otherwise this can prevent the gravitational 
collapse. Fast magnetic reconnection which takes place in the presence of turbulence can induce a 
process of reconnection diffusion. This paper continues our numerical study of this process and its 
implications. In particular, extending our earlier studies of reconnection diffusion in cylindrical clouds 
we consider more realistic clouds with spherical gravitational potentials (from embedded stars) and 
also account for the effects of the gas self-gravity. We demonstrate that within our setup reconnection 
diffusion takes place. We have also derived the conditions under which reconnection diffusion becomes 
efficient enough to make an initially subcritical cloud clump to become supercritical and collapse. Our 
results indicate that the formation of a supercritical core is regulated by a complex interplay between 
gravity, self-gravity, the magnetic field strength and nearly transonic and trans- Alfvenic turbulence. 
In particular, self-gravity helps reconnection diffusion and, as a result, the magnetic field decoupling 
from the collapsing gas gets more efficient compared to the case of an external gravitational field. 
We demonstrate that reconnection diffusion is able to remove magnetic flux from collapsing clumps, 
but only a few of them develop nearly critical or supercritical cores, which is consistent with the 
observations. Their formation is restricted to a range of initial conditions for the clouds as follows: 
thermal to magnetic pressure ratios /3 ~ 1 to 3, turbulent to magnetic energy ratios Eturb/Emag ^ 1-62 
to 2.96, and densities 50 < n < 140 cm~^, when considering ^ 25M0, implying total (gas+stellar) 
masses M^^^ < 12OM0. 

Subject headings: diffusion — ISM: magnetic fields — magnetohydrodynamics (MHD) — star: forma- 
tion — turbulence 



1. INTRODUCTION 
It is well known that star formation occurs in dense 



globule- like regions inside g iant molecular clouds (Blitz 



1993 



]Williams et aD[2QQQD . This is due to theiHow 
temperatures and high densities which help the gravita- 
tional force to overcome the outward forces which act to 
prevent collapse. However, in spite of all observational 
and theoretical efforts, it is not yet fully clear how these 
globules form stars. 

Two key ingredients are known to be present in the 
clouds: magnetic fields and turbulence. On one side, 
the magnetic pressure gradient is an important agent 
working against the gravitational collapse of the gas (see 
Mestel & Spitzer 1956, Mestel 1965, Mouschovias 1991). 
The stability of a cloud supported by the magnetic pres- 
sure may be quantified by the mass-to-magnetic flux ra- 
tio, M/^ N/B^ where M is the cloud mass, B is 
the magnetic field, and N is the column density. This 
ratio in principle defines to what extent a static mag- 
netic field can support a cloud against gravitational col- 
lapse (e.g., Nakano & Nakamura 1978, Crutcher 1999, 
2005a, Crutcher & Troland 2008, Heiles & Crutcher 
2005, Falgarone et al. 2008). Considering an initially 
spherical cloud, the critical value of this ratio implies 
Bcr ^ 2.57 X (7V/10^^ cm"^) /iC, for hydrogen gas. Zee- 



man measurements suggest that the mass-to-magnetic 
flux ratios in molecular cloud cores are around 2.0 times 
larger than the critical value for which gravity an d mag- 
netic support balance each other ( Crutcher||2008 [). 

Turbulence is ubiquitous in interstellar medium r] and 
molecular clouds are part of the big cascade. Turbulence 
rules the structuring and fragmentation of the molecular 
clouds in all scales, and probably also plays an important 
role during the different phases of star forma tion (Mac] 
Low fc Kless"eii|[M4l [McKee fc Ostriker||2007D . TheoE^ 
servation of broad line widths ranging from a lew to more 
than 10 times the sound speed indicate that the turbu- 
lent motions are mainly supersonic to transonic in the 
molecular clouds. Turbulence is also roughly in equipar- 
tition with the magnetic fields in the interstellar medium 
and for this reason it is believed to be trans-Alfvenic 
(Elmegreen & Scalo 2004; Heiles & Troland 2005). In 
fact, the inferred ratio between turbulent and magnetic 
energy for the cold neutral interstellar gas is roughly 
1.3 < Eturb/Emag < 1-9 (Heilcs & Troland 2005). The 
formation of structures appears to be related to com- 
pression caused by the compressional motions associated 
with the turbulent cascade. The relation between com- 
pressible and incompressible parts of the cascade de- 
pends on the sonic and Alfven Machs of the turbulent 
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Sz Lazarian 2001, Chepurnov et al. 2010) and CO lines (see Padoan 
et al. 2006, 2009). 
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flow (Cho & Lazarian 2002, 2003). Even if the turbu- 
lence is globally strong enough for supporting the cloud 
against gravity ([Klessen et al.'200Q|^ Mac Low fc Klessen 
120041 1 Vazquez- Semadeni et al. 20051), compressible mo- 



tions could also cause the collapse at small scales, by gen- 
crating high density regions in bot h magnetized (Heitsch 
eFaI]|2QQT] [Nakamura fcTi|2QQ5l liVlelioli fc de lie 
" ^ Pinof" 



mi 



medium (Klessen et al. 



2006[ [Leao et al.T2009|) and non-ma gnetized 



Elmegreen fc Scalo 2004). 



There are several candidate mechanisms tor injection of 
turbulent energy inside a cloud. Internal sources include 
feedback from low-mass and massive stars. The latter 
in particular, can induce turbulence through intense ion- 
i zing radiation, winds, and supernovae explosions (SNe) 
jMcCray fc Snow'"1979'; 'Melioli fc de Gouveia Dal Ph^ 
[2006; Melioh et al. 2006; Leao et al. 200 ^. External can- 
didates also include SNe shocks (Wada & Norman 2001; 
Elmegreen fc Scalo 20 04; Melioli fc de Gou veia Dal Pino 
2006; Mel ioh et al.|2QQ 6; Leao et al. 2009) andthe galac- 
tic spiral structure ([Kob erts,|1969, ,Bonnell et al. 2006). 
All these processes seem to nave power enough for ex- 
plaining the structuring and kinematics of the ISM, and 
could generate th e observed dispersion of the gas veloci- 
ties (Kornreich fc Scalo 2000). Less energetic sources of 
turbulence include proto-stellar jets and winds, expan- 
sion of HII regions, and magneto-rotational instability, 
but their relative i mportance is still uncertain (|Joung| 
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The exact role that turbulence and magnetic fields play 
on clouds collapse and star formation is still highly de- 
bated. An important question in this debate is that, un- 
der ideal MHD conditions, all the magnetic flux should be 
brought together with the collapsing material and then 
the magnetic field in a proto-star should be several orders 
of magnitude higher than the observed in T-Tauri stars. 
This is often referred as the "magnetic flux problem" 
(see Galli et al. 2006, Johns-Krull 2007, Santos-Lima 
et al. 2012). In order to solve this problem one has to 
consider diffusive mechanisms that can remove magnetic 
flux from a collapsing cloud and allow star formation, at 
least under some circumstances. 

To address this problem of magnetic field diffusion 
both in a partially ionized ISM and in molecular clouds, 
researchers usually invoke the ambipolar diffusion (AD) 
mechanism. AD, which was first discussed in the con- 
text of star formation by Mestel fc Spitzer (1956), has 
been extensively discussed since then (e.g., Spitzer 1968; 
Nakano fc Tademaru 1972; Mouschovias 1976, 1977, 
1979; Nakano fc Nakamura 1978; Shu 1983; Lizano fc 
Shu 1989; Fiedler fc Mouschovias 1992, 1993; Li et al. 
2008; Fatuzzo fc Adams 2002; Zweibel 2002). In princi- 
ple, AD allows magnetic flux to be redistributed during 
the collapse of low ionization regions as the result of the 
differential motion between the ionized and the neutral 
gas. Recent advances in the theory, however, have been 
questioning the efficiency of this diffusion pr ocess in real 
systems. Severa l authors (Shu et al. 2006, Krasnopol 



sky et al. 2010 2011) explored the accretion phase in 



low-mass star formation and concluded that there should 
exist an effective diffusivity about three orders of magni- 
tude larger than the Ohmic diffusivity in order to allow 
an efficient magnetic flux transport to occur. They found 
that AD could work in principle, but only under special 



circumstances, considering specific dust grain sizes (see 
also alternative views in Li et al. 2011). In other words, 
there is still no consensus if AD alone is high enough to 
solve the magnetic flux transport problem in collapsing 
flows. 

Lately, we have been exploring a new mechanism for 
magnetic flux transport whose grounds were originally 
proposed by Lazarian (2005) which is based on the model 
of fast magnetic reconnection in Lazarian fc Vishniac 
(1999). In the presence of magnetic reconnection mag- 
netic fields are not any more frozen irj^and undergo the 
process which was termed reconnection diffusion (see also 
Lazarian 2011; Lazarian et al. 2012b; Santos-Lima et al. 
2010; de Gouveia Dal Pino et al. 2011; 2012; Santos- 
Lima, de Gouveia Dal Pino fc Lazarian 2012, Lazarian, 
Esquivel fc Crutcher 2012). Note, that the problem we 
face in star formation with the magnetic field is not of 
dissipation of magnetic flux, where, indeed, ordinary 
resistivity is necessary, but of magnetic diffusion. The 
ideal non-diffusive MHD condition holds in the absence 
of magnetic reconnection, but when the latter is present, 
changes in the magnetic field topology occur. The model 
by Lazarian fc Vishniac (1999) was confirmed by numeri- 
cal high resolution MHD simulations (Kowal et al. 2009, 
2012) and this provides the justification of our inter- 
preting of our numerical results which do not provide as 
good resolution of current sheets as those in Kowal et al. 
studies. A detailed review of the theoretical foundations 
of this diffusion mechanism by turbulent reconnection is 
presented in the Appendix [A] 

The first numerical study of reconnection diffusion was 
performed in Santos-Lima et al. (2010; see also Lazarian 
et al. 2010; de Gouveia Dal Pino et al. 2011) where we 
explored numerically the role of turbulent reconnection 
on the transport of magnetic field flux from the central, 
denser regions of a molecular cloud to the outer regions. 
This study performed high-resolution three-dimensional 
(3D) MHD simulations of ISM clouds considering a cen- 
tral gravitational field provided by embedded stars. In 
the presence of turbulence, these simulations revealed 
a decrease of the magnetic fiux-to-mass ratio with in- 
creasing density at the center of the gravitational po- 
tential caused by the transport of the magnetic flux to 
the outskirts of the cloud by turbulent reconnection dif- 
fusion. The comparison of these results with systems 
without turbulence and in these cases revealed no change 
in the magnetic-flux-to-mass ratio, as one should expect 
for ideal MHD systems. This was a further evidence that 
in the cases with turbulence the transport is due to re- 
connection diffusion. This effect was observed both when 
starting with equilibrium distributions of gas and mag- 
netic field and when following the evolution of dynami- 
cally unstable (gravitationally collapsing) configurations. 
Thus the process of turbulent magnetic field removal is 
in principle applicable both to quasi-static subcritical 
clouds and collapsing supercritical ones. The increase 
of the gravitational potential, as well as the decrease of 
magnetization of the gas showed an increase of the seg- 
regation of the mass and magnetic flux in the saturated 
final state of those simulations, supporting the notion 

^ The point of the violation of magnetic frozen in condition 
within turbulent plasmas is discussed in detail in Eyink, Lazar- 
ian & Vishniac (2011). 



Turbulent magnetic reconnection 



3 



that turbulent diffusivity relaxes the magnetic field + 
gas system in the gravitational field to its minimal energy 
state. These results were also found to be insensitive to 
the numerical resolution. Tests made with resolutions of 
128^, and 512^ gave essentially the same results as those 
of 256^ thus confirming the robustness of the results and 
the method. 

The study above, however, assumed for simplicity 
clouds with cylindrical gravitational fields. We here 
perform 3D MHD high resolution numerical simulations 
considering more realistic subcritical clouds with cen- 
tral spherical gravitational potentials (due to embedded 
stars) and also account for the effects of self-gravity. 
Our results essentially confirm the results of the previ- 
ous study above, but also reveal the conditions under 
which turbulent reconnection diffusion is efficient enough 
to make an initially subcritical cloud to become super- 
critical and collapse (see also de Gouveia Dal Pino et al. 
2012). 

In Section [2j we describe the setup and the numerical 
methodology that we employed for performing 3D MHD 
simulations of molecular cloud cores or clumps. In Sec- 
tion [3j we present the results of the diffusion of magnetic 
flux by turbulence considering molecular clouds with self- 
gravity and a central spherical gravitational potential. In 
Section [4j we discuss our results and compare them with 
previous works and observations, and in Section [5] draw 
our conclusions. 

2. PHYSICAL SETUP AND NUMERICAL METHODOLOGY 

The astrophysical environment we want to investigate 
in this work consists of an initially subcritical molecular 
cloud clump with a small group of embedded stars sus- 
tained by magnetic fleld and turbulence. Our goal is to 
examine, by means of three-dimensional MHD numeri- 
cal simulations, the conditions under which the trans- 
port of magnetic flux by turbulent reconnection will al- 
low the contraction of the self-gravitating gas to form a 
supercritical core. We also consider a few models with- 
out including self-gravity in order to compare with the 
previous study of Santos-Lima et al. (2010). In most 
of our numerical experiments, the system starts already 
out-of-equilibrium between gravity and the other forces. 
An initially spherical gas clump with a central gravita- 
tional potential that mimics a small group of embedded 
stars is put in the middle of a homogeneous magnetized 
background, then an initial contraction of the gas takes 
place for a short period of time (of the order of the free- 
fall time), after which due to the presence of the mag- 
netic fleld and the turbulence injected, the system evolves 
more smoothly. As an example. Figure 1 shows the log- 
arithmic density map of the central slice of model N2b 
(with no turbulence injected yet) right after an initial fast 
contraction and the injection of turbulence (~ 1.1 Myrs). 
The magnetic fleld conflguration is superimposed to the 
density map. 

We also present one test with the clump starting in 
magneto-hydrostatic equilibrium for comparison with the 
other models. (In this case, we actually started with a 
uniform density cloud, let it relax to the equilibrium and 
then injected turbulence). 

Each system is simulated inside a cubic domain with 
periodic boundaries. Since a typical giant molecular 
cloud has several clumps the use of periodic boundaries 



to model one of these clumps is appropriate. For the sake 
of simplicity, we employ an isothermal equation of state, 
with a single temperature for the whole system, which 
means that it radiates quite efficiently. 




Fig. 1. — Example of a cloud setup in the computational domain 
(with no turbulence injected) right after an initial fast contraction 
at (~ 1.1 Myrs). Logarithmic density map and magnetic field 
lines of the central slice (projected onto the walls of the cubic 
computational domain) as in model N2b of Table 1 with Mpot = 
4O.7M0, initial density no = 90 cm~^, and /3 = 3.0. 



We chose characteristic values for the physical param- 
eters of our models based on observational data (see, e.g. 
Mac Low & Klessen 2004) The computational domain 
has sides L = 3.25pc. For a grid resolution with 256^, 
the corresponding cell length is ^ 0.013pc~ 2600 AU. As 
remarked in Section 1, in the previous numerical study of 
3D MHD turbulent reconnection diffusion in cylindrical 
clouds by Santos-Lima et al. (2010), tests with different 
resolutions between 128^, and 512^ were performed and 
were found to produce essentially the same results, there- 
fore in most of the simulated models here we considered 
a grid resolution of 256^ (see however Fig. ^ below). 

The isothermal sound speed is Cg = 3x lO^m/s, which 
corresponds to a temperature T = 10.9 K (for a mean 
molecular weight ja = 1.0). In order to cover a param- 
eter space as much appropriate as possible to the ob- 
served conditions we considered an initial uniform gas 
density, po, spanning between 1.67 x 10~^^ g cm~^ (or 
no = 10 cm"^) and 1.67 x 10"^^ g cm"^ (or no = 100 
cm~^). We have also considered values of the initial 
ratio between the thermal and the magnetic pressures, 
P = clp/{B'^/S7r) = 3.0, 1.0 and 0.3. This implies initial 
Alfven velocities va ^ 2.4 x lO^cm/s, 4.2 x lO^cm/s, and 
7.8 X lO^cm/s, respectively, and initial uniform magnetic 
field intensities between Bq = 0.35 /iG and 3.4 jaG. 

The initial mass in stars in the cloud clump was 
also varied between 20—60 M©. They produce a spherical 
gravitational potential, ^★(r), given by: 

Mr<w) = -^^ (1) 

^ We note that in this study since we are considering self-gravity, 
the results are not scale independent, therefore, a physical param- 
eter space must be considered 
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where r is the distance to the center of the computa- 
tional box, and 

I max 

1.45pc is a cut-off in the gravita- 
tional force to ensure the symmetry while using period- 
ical boundaries, and is a smoothing radius. Its value 
gives the characteristic length of the clump. We fixed 
= 0.325pc in order to limit the values of the gravita- 
tional force and prevent the system to become Parker- 
Ray leigh- Taylor unstable. 

In order to follow the evolution of the system we solve 
numerically the ideal MHD equations: 



dp 

dt 



dt 



V ■ V V 



V • (pv) = 



1 



(2) 



c^V,+ -(VxB)xB 



gas 



dB 

'dt 



= V X (v X B) 



(4) 



where the independent variables have their usual defini- 
tion. The magnetic field B, satisfies the divergence free 
condition V • B = 0. The potential "^gas is due to the 
gas self-gravity which obeys the Poisson equation: 



(5) 



The source term f of the second equation is the bulk 
force responsible for the injection of turbulence. We 
employ an isotropic, non-helical, solenoidal, delta cor- 
related in time forcing /. It induces harmonic velocity 
fields isotropically distributed in the Fourier space, con- 
centrated around a typical wavelength that defines the 
injection scale linj. The power supply is constant, keep- 
ing the random velocity of the gas Vrms approximately 
constant. In all our models with injection of turbulence, 
we use linj ~ 1.3 pc and Vrms between 2.4 x 10^ and 
5.7 X 10^ cm/s. The injection of turbulent energy starts 
at t = 0.0 Myr and grows slowly until it reaches its max- 
imum at t ~ 10.6 Myrs. 

We employed a modified version of the Cartesian 
Godunov-MH D code o riginally developed by G. Kowal 
(iKowal etal. 12 007; Fal ceta-Goncalves et al. 2008; Santos- 
Lima et al.|20 10), using the HLL solver and Runge-Kutta 
of second order for time integration. The Poisson equa- 
tion is solved by a mu ltigrid based algorithm (Press, 
Teukolsky fc Flannery| p978 )). The code is available 
upon request directly to the authors. 

2.1. Initial Conditions 

Table [T] shows the initial conditions for all the simu- 
lated models both with and without self-gravity. Models 
without self-gravity are labelled " R" , the self-gravitating 
ones starting out of equilibrium are labelled "N" and 
the one starting in magneto-hydrostatic equilibrium is 
labelled "E" model. The "R" models, though scale in- 
dependent, are also described here by their initial condi- 
tions in physical unities for comparison with the other 
models, except for the non-dimensional parameter A 
which measures the strength of the gravitational poten- 
tial A = {GM^)/{L.cl), where G is the Newton gravita- 
tional constant, is the stellar mass potential, L is the 



size of the computational domain (which also gives the 
distance unity in the code), and is the sound speed 
(which is also the velocity unity in the code). The other 
parameters of the Table are the total mass Mtot^ i-e., the 
sum of the mass in stars and the mass in gas, the initial 
cloud gas numerical density no, the initial magnetic field 
^0, the thermal to magnetic pressure ratio /3, the tur- 
bulent velocity '^turb, the Alfvenic velocity va, the initial 
cloud Jeans mass in the presence of magnetic field mj^^, 
and the initial turbulent to magnetic energy density ra- 
tio Efurb/ ^mag foi" the corc < 0.3pc). The turbulent 
and Alfven velocities are given in units of the isothermal 
sound speed Cg = 3 x 10^ cm/s, which was chosen as the 
velocity unity in the code (see below). 

Table |2] shows the dynamical time, tdyn = L/cg. This 
is compared with the free-fall time of the cloud gas which 
is given by. 



(6) 



32 G p 

and the reconnection diffusion time of the magnetic field 



^diff 



V 



(7) 



where the turbulent reconnection diffusivity of the gas is 
(Santos-Lima et al. 2010; Lazarian 2006; 2011; Lazarian 
et al. 2012b): 



if ^turb > va , 



V ^ /inj^turb II ^turb < VA , 

where l^nj = L/kf and '^turb = '^rms- The relations above 
indicate that the ratio ('^turb/'^A)^ is important only in a 
regime of sub- Alfvenic turbulence, i.e. with the Alfvenic 
Mach number Ma < 1. In [7| above, / is the scale where 
the magnetic field is smoothed and is a fraction of the 
size of the computational domain L. In order to compare 
the models, in all simulations with turbulence we consid- 
ered the same 77 (and therefore also the diffusion time). 
Thus for a given /, when changing P we also changed 
Vturb in order to keep 77 constant. This ensured that the 
turbulence decayed at approximately the same time scale 
in all models. 

Also, in order to be able to detect the transport of 
magnetic fiux during the dynamical evolution of the col- 
lapsing clouds, we made the initial turbulent reconnec- 
tion diffusion time smaller than the initial free-fall time 
for all simulated models. We followed the evolution of B 
and p for at least ten time steps, corresponding to ten 
dynamical times, or in total ~ 100 Myrs. For compari- 
son with the turbulent MHD models, we also performed 
simulations with similar initial conditions but without in- 
jecting forced turbulence (these will be hereafter referred 
as "laminar" or non-turbulent models). 

3. RESULTS OF THE NUMERICAL SIMULATIONS 

3.1. Models with self- gravity 

We present here the results for nine models where we 
have included the effects of self-gravity for which the ini- 
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TABLE 1 
Initial conditions for the models 



Model M^{Mq) a no(cm-3) Bo(/iG) /3 vturb{cs) va{cs) mj,s(M0) Eturb/E, 



Rl 


61.1 


0.9 


10 


0.35 


3.0 


0.8 


0.8 


373.7 


2.96 


R2 


40.7 


0.6 


90 


1.06 


3.0 


0.8 


0.8 


314.1 


2.90 


R3 


27.1 


0.4 


100 


1.12 


3.0 


0.8 


0.8 


326.6 


2.89 



Model 


M^Mo) 


Mtot (Mo) 


no (cm 3) 






'^ttirb(Cs) 


Va(Cs) 


rnj,s(M0) 


Eturb/ Emag 


Nl 


61.1 


69.5 


10 


0.35 


3.0 


0.8 


0.8 


373.7 


2.96 


N2a 


40.7 


49.1 


10 


0.35 


3.0 


0.8 


0.8 


448.1 


2.96 


N2b 


40.7 


116.4 


90 


1.06 


3.0 


0.8 


0.8 


314.1 


2.90 


N2c 


40.7 


116.4 


90 


1.85 


1.0 


1.2 


1.4 


379.8 


2.15 


N2d 


40.7 


116.4 


90 


3.37 


0.3 


1.9 


2.6 


639.7 


1.62 


N2e 


40.7 


108.0 


80 


1.74 


1.0 


1.2 


1.4 


389.9 


2.16 


N3 


27.1 


111.2 


100 


1.12 


3.0 


0.8 


0.8 


326.6 


2.89 


N4 


20.4 


104.5 


100 


1.12 


3.0 


0.8 


0.8 


339.6 


2.89 


El 


40.7 


46.6 


90 


1.06 


3.0 


0.8 


0.8 


390.1 


2.90 



TABLE 2 
Characteristic timescales. 



Model 


^dyn(Myrs) 


t//(Myrs) 


tdiffiMyrs) 


Rl 


10.6 


5.7 


1.4 


R2 


10.6 


4.4 


1.4 


R3 


10.6 


4.5 


1.4 


Nl 


10.6 


5.7 


1.4 


N2a 


10.6 


6.7 


1.4 


N2b 


10.6 


4.4 


1.4 


N2c 


10.6 


4.4 


1.4 


N2d 


10.6 


4.4 


1.4 


N2e 


10.6 


4.5 


1.4 


N3 


10.6 


4.5 


1.4 


N4 


10.6 


4.6 


1.4 


El 


10.6 


5.6 


1.4 



tial conditions are listed in Table [T] All models, but 
N2c, N2d, and N2e, have a turbulent velocity which 
is equal to the initial Alfven speed, while these three 
models are sub-Alfvenic, i.e., have a turbulent veloc- 
ity smaller than the Alfven speed. The non-equilibrium 
models Nl and N2a were initialized with the same gas 
density p = 1.67 x 10~^^g/cm ^, but different gravita- 
tional potentials (i.e., different values of the stellar mass 
potential M^). Model N2b has the same initial stellar 
potential as in model N2a, 40. TM©, but a larger gas den- 
sity p = 1.503 X 10~^^g/cm ^ which allows to compare 
the effects of self-gravity. Models N2c and N2d have the 
same initial gas density as in the previous model, N2b, 
but different values of /3 which allow to compare the ef- 
fects of the magnetic field. Model N2e has the same ini- 
tial j3 as in model N2d but a slightly smaller gas density 
p = 1.336 X 10~^^g/cm ^ which allows turbulence to be 
more dominant than in the previous case. The models 
N3 and N4 have the same initial /3 = 3.0 and gas density 
p = 1.67x 10~^^g/cm ^, but different stellar gravitational 
potentials, 27.1M0 and 2O.4M0, respectively. Both have 
smaller stellar gravitational potentials but larger gas den- 
sity (therefore more significant self-gravity) than the pre- 



vious models. Finally, the model starting in equilibrium 
El has conditions which are similar to those of model 
N2b. The remaining initial conditions are the same for 
all models. 

The left panels of Figures [2j [4j and [7| and top- 
left panel of Figure [6j show logarithmic density maps 
of the the central slices of the simulated models after 
100 Myr or ~ten dynamical times. The middle pan- 
els compare the time evolution of the average magnetic 
field-to-density ratio inside a central sphere with radius 
r = 0.3 pc which represents the core of the cloud clump, 
( (5)0.3/ (p)o.3)/(^/p), both for the models with turbu- 
lence (red-dashed line) and without turbulence (black- 
solid line). (The brackets mean averages taken inside the 
cloud radius 0.3 pc.) This ratio is normalized by the av- 
erage value {B/p) taken over all the computational box. 
For the adopted geometry with an initial uniform mag- 
netic field, the ratio above is approximately equal to the 
magnetic flux-to-mass flux ratio within the region con- 
sidered. 

The right panels of these figures and the top-right panel 
of Figure [6] depict profiles of the average magnetic field- 
to-density ratio as a function of the radial distance r from 
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the center of the cloud for the same models, at t = 100 
Myr. 

In Figure |2j we note that after an initial rapid decrease 
caused by the relaxation of the system, the magnetic 
field-to-density ratio inside the cloud cores (middle pan- 
els) remains nearly constant in the laminar cases (i.e., 
without turbulence), while in the cases with turbulence 
there is a clear decrease of it This result indicates 
that there was magnetic fiux transport from the denser, 
more massive central regions to the less dense regions 
outside of the cloud cores. In Figure [2j this effect is par- 
ticularly more pronounced in model N2b (bottom panel) 
which has much larger initial gas density than the other 
models of the Figure, and therefore is under the infiu- 
ence of larger self-gravity. Comparing models Nl and 
N2a which have the same initial gas density, we note 
that the one with larger stellar potential (model Nl in 
top panel) shows a larger decoupling between the mag- 
netic fiux and the mass density. Consistently, the right 
panels of Figure [2] show that the turbulent model with 
larger initial density (for which the effect of self-gravity is 
more important, i.e., model N2b), has a larger decrease 
in the radial profile of the magnetic field-to-density ratio 
in the core region (accompanied by a larger increase in 
the outer regions of the cloud) than the other turbulent 
models, therefore confirming the trend detected in the 
middle panels. 

In order to verify if the initially subcritical clouds have 
developed supercritical cores after the action of the tur- 
bulent reconnection transport, FigurespT] [3l [5]and[8l and 
the bottom panels of Figure [6] compare the radial profiles 
of the integrated magnetic fiux (left panels) , mass (center 
panels), and mass-to-magnetic fiux ratio normalized by 
the critical value (right), at t = 100 Myr for the models 
with and without turbulenc^D 

Figures [4] and [5] show the density maps and similar 
profiles as those of Figure |2] for the models with sub- 
Alfvenic turbulence N2c, N2d, and N2e. Models N2c 
and N2d have the same initial gas density and stellar 
mass as in model N2b, but smaller initial (3 and turbu- 
lent to Alfven velocity ratio. This naturally makes the 
decoupling of the magnetic fiux from the dense gas re- 
gions by turbulent reconnection diffusion more difficult. 
Both models N2c and N2d, for which p = 1.0, and 0.3, 
respectively, have no significant magnetic fiux transport 
outward when compared to model N2b in Figure [2] for 
which (3 = 3.0, or to their laminar counterparts. How- 
ever, when the initial mass density of model N2c is de- 
creased from 90 to 80 cm~^ (keeping P = 1.0), as in 
model N2e, the transport of magnetic fiux by turbulence 
is slightly enhanced (though still much less than in model 
N2b) and the initially subcritical cloud develops a nearly 
critical core after 100 Myr (see Figure [5|. Therefore, this 
result indicates that the increase of gas density and to- 
tal gravity will not always result in enhancement of the 
fiux transport by the turbulence as we have seen in the 
models of Figure [2] When turbulence is sub-Alfvenic 

^ We note that the oscillations observed in these plots (which are 
slightly stronger in the laminar models) are acoustic oscillations of 
the cloud due to the fact that the virialization time of these systems 
is larger than the simulated period. 

^ The critical value of the mass-to-magnetic flux ratio for a disk 
shaped system is given by: {27tG-^^'^)~-^ 



(and thus the magnetic field is strong for the level of tur- 
bulence applied), a smaller gravitational potential will 
delay the collapse and thus may give time for the tur- 
bulent reconnection diffusion to transport outward part 
of the magnetic flux. This explains why the decrease in 
density (and this of self-gravity) from no = 90 cm~^ in 
model N2c to no = 80 cm~^ in model N2e, results in the 
build up of a nearly critical core in the last case. 

Figure |6] shows the density map and proflles for model 
El. This model which starts in magneto- hydrostatic 
equilibrium has the same initial f3 and density po in its 
central region as the non-equilibrium model N2b, but the 
equilibrium condition makes its density proflle stratifled 
and thus the total mass inside the cloud is smaller than 
in N2b model (see Table [l]). As in the turbulent N2b 
model. El also suffers outward magnetic flux transport, 
as indicated by the decreasing magnetic fleld-to-density 
ratio in the core region (top- middle panel). However, the 
much smaller total mass of El prevents it to collapse and 
become a supercritical core, as indicated by the mass-to- 
magnetic flux diagrams of Figure [g] (bottom-right panels, 
red-dashed line). 

Figures W\ and [8] show the results for models N3 and 
N4 which have the same initial (3 = 3.0 and gas den- 
sity no = 100 cm~^ (which are comparable with model 
N2b for which the initial gas density is only a factor 
0.1 smaller), but different stellar gravitational potentials, 
27.IM0 and 2O.4M0, respectively (which are smaller 
than that of model N2b). 

In both turbulent models we note a smaller decrease 
with time of the magnetic fleld-to-density ratio in the 
cloud core relative to model N2b (see middle panels of 
Figure I?]), specially for the model with smaller stellar 
potential (N4). The inspection of the radial proflle of 
the mass-to-magnetic flux ratio at t = 100 Myr for this 
model in Figure |8] shows that in fact, there is no trans- 
port of magnetic flux in this case. This occurs because of 
the combination of two effects. The larger magnetic fleld 
of model N4 relative to model N2b (see Table fl]) makes 
it more diflicult for the turbulence to decouple the mag- 
netic flux from the denser material. At the same time, 
the smaller stellar gravitational well slows down the in- 
fall of matter to the center in model N4 (see Table [2|, 
while the action of the turbulence helps to spread the gas 
against gravity (bottom left panel of Figure [8|. Model 
N3 on the other hand, although with the same strength 
of magnetic fleld as model N4, it has a stellar potential 
large enough to push the core matter inwards while tur- 
bulence decouples the magnetic flux transporting it to 
the less dense regions, allowing the formation of a su- 
^rcritical core (as we see in the upper panels of Figure 

3.2. Effects of Resolution on the Results 

As remarked, all the numerical simulations presented 
above have been performed with a resolution 256^. In 
order to ensure that the results above are not being af- 
fected by numerical effects, we also ran one of the models 
(model N2e, Table [T]) with increased and decreased reso- 
lutions. Fig. [9] compares the same quantity presented in 
Figure |4] (middle panel) for this model with resolutions 
128^, 256^, and 512^. We do not observe signiflcant dif- 
ference between them. Thus, we can expect that the 



Turbulent magnetic reconnection 



7 




250 1 



1501 



1501 




0,0 



0,5 



1.0 

(pc] 



1,5 



2,0 



40 60 

t (MvO 

Fig. 2. — Model Nl with Mpot = GLIM© (top), models N2a (center) and N2b (bottom) with Mpot = 40. TM©. Top and middle models 
have initial densities no = 10 cm~^, and the bottom model N2b has no = 90 cm"*^. All models have (3 = 3.0 Left panels show logarithmic 
density maps of the central slices of the turbulent cloud models at t = 100 Myrs. Middle panels show the temporal evolution of the 
average magnetic field-to-density ratio at the cloud core region of radius rc = 0.3 pc normalized by the average value over the entire cloud, 
({B)o_s/ {p)o_s)/iB/ p), for the turbulent (red-dashed lines) and the laminar (black continuous lines) models. Right panels show the radial 
profile of the average magnetic-to-density ratio at t = 100 Myrs {B) / {p) for these models. 

results presented for the models with resolution of 256^ 
are robust. 



3.3. Models with no self- gravity 

We have not discussed yet the specific effects that the 
inclusion of self-gravity produces on the turbulent core 
collapse. In order to due that, we have to compare the 
self-gravitating models analyzed in the previous section 
with counterparts without self-gravity. These models la- 
belled "R" are also listed in Table [H Model Rl has 
similar initial conditions to those of the self-gravitating 
model Nl of Table 1, with a gravitational potential pa- 
rameter A = 0.9 which is equivalent to a stellar mass of 
61.7 Mq. Models R2 and R3, in turn, have similar initial 
conditions to the self-gravitating models N2b and N3 of 
Table 1, with 40.7 Mq ^ A = 0.6, and = 27.1 

Mq A = 0.4, respectively. 

The logarithmic density maps at 100 Myrs as well as 



all the related radial and time evolution profiles for these 
models are presented in Figures [lO] and [TT] 

The comparison of the models without self-gravity R2 
and R3 with the self-gravitating models N2b and N3, 
respectively (bottom panels in Figures |2] and [3] and up- 
per panels in Figures [t] and [8| reveals the importance 
of the inclusion of self-gravity, particularly in model R2 
which does not develop a supercritical core contrary to its 
self-gravitating counterpart (model N2b) . On the other 
hand, model Rl (top panels in Figures 10 and 11) is not 
much affected by the elimination of sell-gravity, as we 
can see when comparing with the self-gravitating model 
Nl (top panels of Figures [5] and [3| . In this case the ex- 
ternal gravitational potential due to the embedded stars 
is already very high and dominates the cloud collapse 
making self-gravity non-negligible only in a very small 
radius close to the center. 
Model Rl is also comparable to the model D2 of 
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Fig. 3. — Radial profiles at t = 100 Myr for the magnetic fiux <l> (left panels), the mass M (center panels), and for the magnetic fiux-to- 
mass ratio normalized by the critical value (M / ^) / (M / ^) crit (right panels) for the models of Figure |2] Top panels: model Nl (no = 10 
cm~^, and = 61. IM©); center panels: model N2a (no = 90 cm~^, and Mpot = 40. TM©); bottom panels: model N2b (no = 10 cm""^). 
Red-dashed lines are for turbulent models and black continuous lines are for the laminar models. 

Santos-Lima et al. (2010). Both have the same initial 
conditions and no self-gravity, but differ in the geometry, 
as the clouds in Santos-Lima et al. have initial cylindrical 
gravitational potential rather than spheric. The presence 
of a spherical gravitational field yields a smaller efficiency 
in the magnetic flux transport, as we can see when com- 
paring the top- middle diagram of Figure |2] with top-right 
panel of Figure 11 of Santos-Lima et al. (2010). This is 
because, while in the case of a spherical potential all mat- 
ter is pushed to a single central point, in the case of a 
cylindrical potential, gravity pushes the collapsing mate- 
rial to the central axis all along the cylinder and thus it 
is more effective to help the decoupling between the col- 
lapsing gas and the magnetic flux driven by the turbulent 
reconnection diffusion. 

4. DISCUSSION 

We have performed 3D MHD simulations of the evolu- 
tion of self-gravitating cloud clumps with a central grav- 
itational potential with spherical symmetry and embed- 



ded in an initially uniform magnetic field. For compar- 
ison we have also considered a few models without self- 
gravity. The simulations were started with the cloud 
either in magneto-hydrostatic equilibrium (see model El 
of Table 1) or out of it (all the other models of Table 
1). We injected non- helical turbulence in each system 
with an rms velocity (vrms) comparable both to the ini- 
tial Alfven speed and to the isothermal sound speed (see 
Table [T]) and then let it to evolve. All the tested models 
have an initial turbulent to magnetic energy ratio > 1, 
which is compatible with estimates from observations of 
the cold neutral interstellar medium (Heiles & Troland 
2005). We have imposed to all systems the same turbu- 
lent time decay (eqs. ^ in order to compare the effects of 
the turbulent reconnection diffusion among them. Also 
for comparison, we evolved the same set of models with 
no injection of turbulence. 

Most of the turbulent models here investigated evi- 
dence the transport of magnetic flux, i.e., the decoupling 
of the magnetic field from the denser, inner regions of 
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Fig. 4. — The same as in Figure |2] for models N2c, N2d, and N2e. Top and Bottom: models N2c and N2e with /3 = 1.0 Center: model 
N2d with /3 = 0.3. Top and middle models have initial densities no = 90 cm~^, and the bottom model has no = 80 cm""^. All models have 



the cloud clump due to the presence of MHD turbu- 
lence. The exceptions are model R3 without self- gravity 
(as demonstrated by Figures 10 and 11) and the self- 
gravitating models N2c, N2d, and N4 (see Figures [4] and 
[sj and[7|and[8j respectively). 

4.1. Comparison with Santos-Lima et al. (2010) results 

Previous numerical studies of cylindrical cloud systems 
by Santos-Lima et al. (2010; see also Lazarian 2011; de 
Gouveia Dal Pino et al. 2011) had already evidenced the 
importance of the effects of turbulent magnetic reconnec- 
tion to remove the magnetic flux excess from collapsing 
systems, as originally suggested by Lazarian (2005) based 
on the fact that magnetic reconnection must be fast in 
weakly MHD turbulent environments (Lazarian & Vish- 
niac 1999; Kowal et al. 2009; see also Appendix A). The 
present numerical study, which considers more realistic 
self-gravitating spherical clouds, confirms the results of 
this previous analysis, i.e., it shows that the presence of 
turbulent reconnection diffusion is able to remove mag- 
netic flux from the denser central regions to the edges 



of the cloud, therefore facilitating the gravitational col- 
lapse in most of the tested models, without considering 
ambipolar diffusion effects. This is assured by the mea- 
sured magnetic-flux-to-mass ratio in the core regions of 
the simulated clouds, which is quantified in our study 
by the average magnetic field-to-density ratio along the 
magnetic field lines (see middle diagrams of Figures |2j 
[4j [6]-top, j7| and[lO|. This ratio decreases with time m 
most of the turbulent models, while it remains constant, 
on average, in the non-turbulent counterpart models. 

As in Santos-Lima et al. (2010), we also find that, in 
general, an increase in the stellar gravitational potential 
(see e.g., models Nl and N2a of Figure [2|, as well as a 
decrease in the initial magnetic field strength (or an in- 
crease in B\ see models N2b of Figure [2] and N2c and N2d 
of Figure p| favours the turbulent reconnection transport 
of the magnetic flux and its decoupling from the dense 
collapsing gas. 

In the cases when the turbulence is sub-Alfvenic (i.e., 
"^turhhA < 1, as in models N2c, N2d, and N2e), the flux 
transport by turbulent reconnection diffusion is more dif- 
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ficult, as one should expect, since the turbulence at large 
scales is weak. In fact, for a cloud with an initially too 
large magnetic field, sub-Alfvenic turbulence will not be 
able to transport the magnetic fiux to outside and the 
cloud may fail to build up a supercritical core, i.e., a 
core with a mass-to-magnetic fiux ratio above the crit- 
ical value that is required for gravity to overcome the 
magnetic forces (see for instance, models N2c and N2d in 
Figures Hand [5] for which /3 = 0.3 and 1.0, and Vturb/'^A = 
0.7 and (J.9, respectively). 

Nonetheless, even sub-Alfvenic regimes of turbulence 
may allow the formation of critical cores. This was the 
case of model N2e, which has the same initial conditions 
of model N2c, except for a smaller total gravitational 
potential. The latter caused a delay of the gas collapse 
that gave time for the sub-Alfvenic turbulent reconnec- 
tion, which becomes stronger at smaller scales (Lazarian 
2006), to transport outward part of the magnetic fiux 
and allow the formation of a critical core (see Figure [5|. 

This result can be understood in terms of the recon- 
nect ion diffusion coefficient in sub-Alfvenic regimes. Al- 



though weak at large scales, as the turbulence cascades 
the strength of the interactions increases and at a scale 
/ linj{vturb/vA)'^ it bccomcs strougcr and therefore, 
more efficient to help with magnetic ffux transport. As 
stressed by Lazarian (2006, 2011) and also tested nu- 
merically by Santos-Lima et al. (2010), the diffusivity 
in this regime is expected to be given approximately by 
hnj'^turhiyturh / (^^^ ^Q- ^1 ^hus Smaller than the hy- 
drodynamic turbulent diffusivity by a factor iyturh/vAY ^ 
with the eddies of the strong turbulence playing a critical 
role. In the next section we discuss the validity of this 
estimate for the turbulent reconnection diffusion coeffi- 
cient. 

4.2. Comparison of magnetic reconnection diffusivity 
and resistivity effects 

As in Santos-Lima et al. (2010), we can evaluate 
the effective turbulent reconnection diffusion coefficient 
of our simulated models by comparing them with non- 
turbulent resistive models with enhanced Ohmic resis- 
tivity. For instance, considering the same initial condi- 
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tions as those of the sub-Alfvenic models N2c, N2d and 
N2e, we performed several simulations of resistive non- 
turbulent models considering different values of Ohmic 



resistivity. Table 3 presents a set of resistive models 
N2cr, N2dr, and N2er whose initial conditions are the 
same as those of N2c, N2d and N2e, respectively. Figure 
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Fig. 8. The same as in Figure [s] for models N3 and N4 of Figure [t] Top panels: model N3 (no = 100 cm-^, and = 27.IM0); 
bottom panels: model N4 (no = 100 cm~3, and = 2O.4M0). 

tions although further developments are still required in 
order to obtain a more precise values of the reconnection 
diffusion coefficients (see Lazarian 2011 and references 
therein). 

TABLE 3 

Parameters for resistive models with same initial 

CONDITIONS OF THE SUB-AlFVENIC MODELS N2C, N2D AND N2E 




Fig. 9. — Comparison between different resolutions for the model 
N2e (Table [ij. 



12 compares model N2c with two resistive models N2crl 
and N2cr2 with same initial conditions and enhanced 
Ohmic resistivities r]ohm = 0.005 c.u. and r]ohm = 0.001 
C.U., respectively. We see that the resistive model with 
Vohm = 0.001 is the one which best reproduces model 
N2c, thus we can conclude that the effective turbulent 
resistivity of model N2c is rjturb ^ 0.001 c.u. A similar 
procedure for models N2d and N2e allowed us to estimate 
also their turbulent diffusion coefficients, which are pre- 
sented in Table [4] This Table also shows the correspond- 
ing upper limit estimates for the ratio rjturb/hnj'^turb (and 
the Alfvenic Mach number) for these models. According 
to eq. |8l we should expect this ratio to be smaller than 
1 for sub-Alfvenic turbulence. Therefore, the obtained 
values are roughly consistent with the theoretical predic- 



Model 


n (cm ^) 




VOhm (c.u.) 


N2crl 


90.0 


1.0 


0.001 


N2cr2 


90.0 


1.0 


0.005 


N2cr3 


90.0 


1.0 


0.01 


N2cr4 


90.0 


1.0 


0.5 


N2cr5 


90.0 


1.0 


0.1 


N2drl 


80.0 


1.0 


0.001 


N2dr2 


80.0 


1.0 


0.005 


N2dr3 


80.0 


1.0 


0.01 


N2erl 


90.0 


0.3 


0.0005 


N2er2 


90.0 


0.3 


0.001 


N2er3 


90.0 


0.3 


0.005 



TABLE 4 

The diffusivity coefficient for the sub- Alfvenic models 



Model rjturb (cu.) 


Vturb/ "^turbhnj 


(Vturb/VA)^ 


N2c < 0.001 


< 0.0024 


~ 0.229 


N2d < 0.001 


< 0.0016 


~ 0.096 


N2e < 0.001 


< 0.0020 


~ 0.315 
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Fig. 10. Model Rl with A = 0.9 (Mpot = 61. IM©) (top), model R2 (center) with A = 0.6 {Mpot = 40. TM©), and model R3 (bottom) 
with A = 0.4 (Mpot = 27.IM0). The models have initial densities po = 1.0, and /5 = 3.0. Left panels show logarithmic density maps 
in the xz-plane of the central slices of the turbulent cloud models ait t = 100 Myrs. Middle panels show the temporal evolution of the 
average magnetic field-to-density ratio at the cloud core region of radius Tc = 0.3 pc normalized by the average value over the entire cloud, 
{{B)o.3/ {p)o.3)/{B/ p), for the turbulent (red-dashed lines) and the laminar (black continuous lines) models. Right panels show the radial 
profile of the average magnetic-to-density ratio at t = 100 Myrs (B) / (p) for these models. 



4.3. Cylindrical versus spherical gravitational potentials 

Another result of particular importance here is the fact 
that in the presence of a more realistic spherical gravi- 
tational field in the cloud, the magnetic flux transport is 
less efficient than in the presence of a cylindrical field (as 
in Santos-Lima et al. 2010). This was already expected, 
since in a spherical potential all gas is pushed to a sin- 
gle central point, while in the cylindrical field the gas 
is pushed to the central axis along the cylinder making 
the decoupling of the gas from the magnetic flux, which 
is driven by the turbulent reconnection diffusion, more 
effective. Therefore, we may conclude that the results of 
the previous study based on cylindrical turbulent clouds 
by Santos-Lima et al. (2010) have overestimated the flux 
transport by reconnection diffusion. For instance, for 
model Rl with spherical potential (top panels of Figure 
|2| the normalized magnetic flux-to-density ratio has de- 



creased by a factor 0.11 after 8 dynamical times, while 
the same model with cylindrical field in Santos-Lima et 
al. (2010) (model D2) has decreased by a larger amount, 
0.65, at the same time interval. 

4.4. Effects of self- gravity upon the magnetic flux 
transport by reconnection diffusion 

The comparison of models without self-gravity with 
self-gravitating models have revealed that self-gravity 
can also significantly help the decoupling between gas 
and magnetic flux due to reconnection diffusion, partic- 
ularly in the late stages of the cloud collapse. A critical 
example is the self-gravitating model N2b (Figures |2] and 
[3| in which reconnection diffusion causes the build up of 
a supercritical core, while its counterpart without self- 
gravity, model R2 (Figures p^O]), is unable to develop a 
supercritical core. 

An increase of the self-gravity (which is provided by an 
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Fig. 11. — Radial profiles at t = 100 Myr for the magnetic fiux <l> (left panels), the mass M (center panels), and the magnetic fiux-to-mass 
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= 6I.IM0); center panels: model R2 (A = 0.6 , i.e. = 41. TM©); bottom panels: mo^ R3 (A = 0.4 , i.e. = 27.IM0). 

Red-dashed lines are for turbulent models and black continuous lines are for the laminar models. 



increase in the initial gas density of the cloud) improves 
the turbulent transport of the magnetic flux. If the gas 
density in the cloud is large enough (no > 50 cm~^), 
its effect seems to be more important than that of the 
stellar gravitational potential (for ~ 41 Mq) to help 
the decoupling between the gas and the niagnetic field 
(see models N2b and N3 in Figures [2] and [sTand ItI and 
pi. However, for a given strength 01 burbulence, 11 one 
mcreases the initial density or the stellar gravitational 
field indefinitely, then eventually the total gravitational 
potential will become so high that it will neutralize the 
ability of the reconnection diffusion to decouple the mag- 
netic flux from the dense gas and thus, no efficient trans- 
port of magnetic flux will occur to outside of the cloud 
core, as most of it will be dragged by the infalling gas. 
We have seen this effect, for instance, when increasing 
the initial cloud density in model N2e (which has initial 
/3 1.0, - 41 M0, and no = 80 cm"^) to no = 90 
cm~^ in model N2c (see Figures [4]and[5|. While the 
first model evidences some magnetic flux transport and 



develops a nearly critical core, the second one fails com- 
pletely. Similarly, for model N2b (which has P = 3.0, 
^ 41 M0, and no = 90 cm~^), transonic, trans- 
Alfvenic turbulence allows the formation of a supercriti- 
cal core, as indicated in Figures [2] and [S] However, if one 
increases its initial gas density to 100 cm~^, the total 
gravitational potential becomes so large and makes the 
collapse so fast that the turbulent magnetic reconnection 
becomes ineffective to decouple the magnetic flux from 
the dense core. On the other hand, if one reduces the 
stellar mass of model N2b to ~ 27 Mq and increases 
the gas density to 100 cm~^ as in model N3 (Figures [t] 
andlsl, so that the total mass is nearly the same as in 
model N2b, then some flux transport is evidenced and 
a supercritical core develops, but a further reduction of 
the total mass (as in N4) again prevents the formation 
of a critical core, because in this case the infall becomes 
so slow that the turbulence actually helps to spread out 
the core material. 
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Fig. 12. — Temporal evolution of the average magnetic field-to- 
density ratio at the cloud core region within a radius Tc = 0.3 
pc normalized by the average value over the entire cloud, for the 
turbulent N2e model (black, solid line), the resistive N2erl model 
with riohm — 0.001 c.u. (blue, dot-dashed line); and the resitive 
N2er2 model with rjohm = 0.005 c.u. (red, dashed line). We 
note that the resistive model N2erl is comparable to the turbulent 
model N2e. 

4.5. Effects of the cloud initial conditions 

All the results above were found for cloud clumps which 
had initial uniform density and were out of magneto- 
hydrostatic equilibrium when turbulence was injected. 
We have also tested a model starting in magneto- 
hydrostatic equilibrium having a stratified density, with 
the central density and the remaining initial conditions 
as in our reference model (N2b) (see model El in Fig- 
ure [6|. As N2b, model El also undergoes an efficient 
outward magnetic ffux transport due to turbulent recon- 
nection. However, the much smaller total mass of El due 
to the cloud stratification (with a fraction 0.43 of the to- 
tal mass of N2b model) prevents its collapse to form a 
supercritical core within the evolved time interval. 

The results above indicate that the formation of a su- 
percritical core is regulated by a complex interplay be- 
tween gravity, self-gravity, the magnetic field strength 
and nearly transonic, trans- Alfvenic turbulence. Al- 
though we have found that turbulent reconnection diffu- 
sion is very efficient to remove magnetic ffux from most of 
the collapsing core clump models tested here, only a few 
were succeeded to develop nearly critical or supercritical 
cores (see models N2b, N2e, and N3) which may be able 
to collapse and form stars. In other words, for the cloud 
conditions investigated here, the formation of supercriti- 
cal cores is restricted to a limited range of parameters, as 
one actually should expect from observations that pre- 
dict a low efficiency of star formation (see e.g. Mac Low 
& Klessen 2004; Leao et al. 2009; Vazquez- Semadeni et 
al. 2011). 

To summarize, our results suggest that turbulent re- 
connection ffux transport will allow initially sub-critical 
clouds to become nearly critical or supercritical for cloud 
clumps with initial values of /3 ~ 1 to 3, cloud densities 
50 < no < 100 cm~^ when considering stellar masses 
~ 4IM0, and densities 100 < no < 140 cm"^ when 
considering stellar masses ~ 27M0, implying total 
cloud clump masses Mtot ^ 12OM0. For smaller densi- 
ties the clouds are fragmented by the injected turbulent 



power and no core is build up|^ For higher densities, the 
effects of self-gravity are so strong that the core collapses, 
dragging most of the magnetic ffeld, so that no signiff- 
cant magnetic ffux is detected. Of the 9 self-gravitating 
models here investigated, 4 formed marginally critical or 
supercritical cores (Nl, N2b, N2e and N3) and 2 subcrit- 
ical cores (N2a and El), all of which evidenced turbulent 
magnetic flux transport. The 3 remaining models (N2c, 
N2d and N4) did not evidence any magnetic flux trans- 
port by reconnection diffusion (either because the initial 
cloud had too strong magnetic ffeld or too strong turbu- 
lence, as described in Section 3). Table [s] lists the ffnal 
conditions of the built up cores all simulated models. 

4.6. Comparison of our results with observations 

The direct approach to resolve the uncertainty con- 
cerning the process by which stars form is to measure 
the magnetic ffeld strengths in molecular clouds in or- 
der to see whether they are weak enough (i.e., the cloud 
is supercritical) or strong (i.e., the cloud is subcriti- 
cal). Recent Zeeman measurements of dark cloud cores 
and envelopes by Crutcher et al. (2009, 2010; see also 
Troland & Crutcher 2008) have raised new challenges 
to the current theory of magnetic ffux transport by am- 
bipolar diffusion. They obtained the mass-to-magnetic 
ffux ratio between the cloud core and the envelope, 
R = {Mc/^c)/{Me/^e)^ as well as the mass-to-ffux ratio 
between the core and the entire (core+envelope) cloud, 
= {Mc/^c)/{Mc^e/^c+e), and found that these ratios 
are less than unity for four observed clouds (Bl, B217-2, 
L1448CO and L1544). The ambipolar diffusion requires 
the mass-to-magnetic ffux ratio in the core to be larger 
than in the entire cloud, resulting R' > 1^ which is in 
contradiction with the results found by Crutcher et al. 
Lazarian et al. (2012a) provided the explanation of the 
Crutcher et al. (2009) results that is based on the con- 
cept of reconnection diffusion. They, however, did not 
do numerical simulations and in their calculations relied 
on the phenomenological description of reconnection dif- 
fusion. 

Although the conditions of the cloud core-envelopes 
investigated by Crutcher et al. are somewhat distinct 
from those obtained from the numerical simulations in 
the present study, we can make at least qualitative com- 
parisons with their results. For instance, although denser 
and more magnetized, their cores are all subcritical or 
only marginally critical (i.e., they have mass-to- magnetic 
ffux ratios relative to the critical value in the range 
l^crit = 0.45 — 1.15). The cores built up in our mod- 
els which evidenced turbulent ffux transport have ffnal 
average mass-to-magnetic ffux ratios jj^crit = 0.15 — 5.25. 
Most of these cores have ratios R and R' which are con- 
sistent with the inferred ones by Crutcher et al. (2009), 
with the exception of the cores built up in models Nl 
and Nb2 (see Table 3, where these ratios are presented 
for the ffnal snapshots of all our models). Indeed, obser- 
vationally, it is known that the level of turbulence drops 
in the cores. We have seen that reconnection diffusion 
slows down with the decrease of the turbulent velocity 
(Santos-Lima et al. 2010; de Gouveia Dal Pino et al. 

^ It is possible that this effect has been a little overestimated 
here due to the employment of an isothermal equation of state for 
the clouds. 
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2011; Lazarian 2011), so that as discussed in Lazarian 
et al. (2012a) we should expect a slower transport of 
magnetic flux to outside of the core as compared to the 
envelope. This behaviour is actually detected in the cores 
formed in the present study. 

In the particular case of model N2b, the computed ra- 
tios R and R' are smaller than unity until the core be- 
comes highly supercritical, from this point on R and R^ 
become larger than unity due to the runaway increase of 
the mass-to-flux ratio of the collapsing core (maximum 
l^crit.c ^ 40.7, as shown in Table [5 and Figure l3|. 

Figures 13 and 14 illustrate the transport oiThe mag- 
netic field lines from the collapsing core to the surround- 
ing envelope for model N2b. While Figure 13 compares 
this turbulent model with its non-turbulent counterpart 
at the final snapshot (t = 100 Myr), Figure [m] depicts the 
time evolution of this turbulent model untiTthe start of 
the core collapse around 90 Myr when then, the portion 
of the magnetic flux that was not diffusively transported 
to the surrounding envelope is advected to the center by 
the collapsing gas causing a large increase of the mag- 
netic field intensity in this region. 

In the case of model Nl, the ratios R and R are also 
smaller than 1 for several time-steps. After 68 Myrs, the 
very rapid increase of the central density caused by the 
high central gravitational potential in this model, makes 
R and R to increase to values larger than 1 (the same 
situation applies to the non-self- gravitating model Rl). 

Therefore, the three cores end up with R and R > 1 
because they are already collapsing to form proto-stars, 
N2b due to the dominance of self- gravity and Nl (and 
Rl) due to the strong central potential. The other su- 
percritical cores (N2e and N3) still have values of R and 
< 1 at the final snapshot, therefore comparable to 
those obtained for the observed cores by Crutcher et al. 

We should note that the inferred values for R and R' 
from the observations are subject to significant uncer- 
tainties due to measurement limitations and simplified 
assumptions (see Mouschovias & Tassis 2010). There- 
fore, while encouraging, the results from the comparisons 
above should be viewed with caution. Besides, they call 
for further and more precise observational estimates of 
these ratios (see also Lazarian et al. 2012a). 

5. CONCLUDING REMARKS 

The mechanism here discussed of reconnection diffu- 
sion may present the last missing piece for constructing 
the new paradigm of star formation where turbulence and 
turbulent feedback play a central role (see also Lazarian 
2011). The present numerical study together with the 
earlier one by Santos-Lima et al. (2010) have investi- 
gated this new mechanism focussing on the early stages 
of star formation. Recently, Santos-Lima et al (2012a, 
2012b) have also investigated this mechanism in the late 



stages, during the formation of the protostellar disks. 
These authors have shown, by means of 3D MHD sim- 
ulations, that the diffusivity arising from magnetic re- 
connection in the presence of turbulence is also able to 
transport magnetic flux to the outskirts of the (cloud 
core) disk progenitor at time scales compatible with the 
core collapse. In just a few lO^yr, a rotationally sup- 
ported disk forms around the protostar with a nearly 
Keplerian profile as required by the observations. Since 
MHD turbulence is expected to be present also in pro- 
tostellar disks, this is a natural and fast mechanism for 
removing magnetic flux excess and allowing the forma- 
tion of these disks. Therefore, taken together these stud- 
ies have tested this new paradigm demonstrating that 
magnetic flux removal from molecular clouds by recon- 
nection diffusion is very efficient and calls for reconsider- 
ing the relative role played by ambipolar diffusion (and 
other transport mechanisms) in the processes of star and 
planet formation. 

Finally, we should remark that in this study, we have 
focused on the evolution of isothermal self-gravitating 
clouds with embedded stars which provided an external 
gravitational potential. As a matter of fact, some of the 
mentioned observed dark cloud cores contain stars em- 
bedded in them (e.g., Bl and L1448CO clouds; Bachiller, 
Menten, & del Rio- Alvarez 1990; Volgenau et al 2006). 
In forthcoming work, we will explore the effects of the 
magnetic flux transport by turbulent reconnection in the 
evolution of initially starless clouds in order to assess 
the effects of self- gravity only upon the transport, with- 
out considering an external field. Also, the isothermal 
approximation assumed here mimics the effects of an ef- 
ficient radiative cooling of the gas. However, in more 
realistic cases, a detailed treatment of non-equilibrium 
radiative cooling in the clouds (e.g., Melioli et al. 2005) 
is required, particularly in the late stages of the core for- 
mation. The effects of non-equilibrium radiative cooling 
will be also considered in these forthcoming studies. 
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APPENDIX 



THEORETICAL GROUNDS OF MAGNETIC FLUX TRANSPORT BY TURBULENT RECONNECTION 

The magnetic diffusion mechanism that we address here is a process deeply rooted in the microphysics behavior of the 
magnetic fields in highly conductive flows. Textbooks characterize these flows by the Lundquist number S = LxVa/ti^ 
where t] is the Ohmic diffusivity, is a typical scale of the system and va is the Alfven velocity. For astrophysical 
systems is in general very large and therefore, S ^ 1 which makes magnetic diffusion negligible. However, one 
may ask: does the magnetic field remain absolutely frozen-in within highly ionized astrophysical fluids? The answer 
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Fig. 13. — Magnetic field line distribution and the equatorial cut of the logarithmic density distribution for the collapsing cloud model 
N2b. The right panel shows the turbulent model while the left panel shows its laminar counterpart for comparison. The time depicted is 
100 Myrs for both models. The central core has a density 50 cm~^. 100 iqually spaced magnetic field lines are plotted within a radius of 
~ 1.3 pc. The color scale indicates the intensity of the magnetic field whith red representing the maximum and blue the minimum values. 



TABLE 5 

Final quantities obtained for the built up cores and envelopes. Central density Uc and magnetic field in z-direction Be, 

TURBULENT-TO-MAGNETIC ENERGY RATIO Eturb/Emag, MASS Mc, MAGNETIC FLUX ^c, AND MAXIMUM MASS-TO-MAGNETIC FLUX RATIO 
RELATIVE TO THE CRITICAL VALUE /J,crit,c FOR THE BUILT UP CLOUD CORES (Tc < 0.3 PC). MASS Me AND MAGNETIC FLUX FOR THE 
FORMED CLOUD ENVELOPES; MASS-TO-MAGNETIC FLUX RATIOS BETWEEN THE CLOUD CORE AND THE ENVELOPE, R, MASS-TO-FLUX RATIOS 
BETWEEN THE CLOUD CORE AND THE ENTIRE CLOUD (CORE+ENVELOPE) CLOUD R , AT THE FINAL TIME STEP {t ~ 100 MYRS). MODELS Nl, 

N2b, N2e, N3 and R1 developed supercritical or marginally critical cores. All models but N2c, N2d, and N4, evidenced 

magnetic flux transport by turbulent reconnection diffusion. 



Model 


nc(cm ^) 




Efurb/ Emag 


Me (Mo) 


^c(103iG.cm2) 


f^crit,c 


Me (Mo) 


<l>e(103iG.cm2) 


R 


R 


Rl 


3.8 X 10^ 


1.92 


4.83 


0.7 


0.20 


1.32 


7.7 


3.54 


1.52 


1.47 


R2 


4.2 X 10^ 


4.04 


1.68 


2.0 


0.62 


0.96 


73.7 


10.6 


0.46 


0.48 


R3 


994.7 


3.32 


3.55 


1.1 


0.43 


0.73 


83.0 


11.2 


0.33 


0.34 


Nl 


3.5 X 10^ 


2.60 


4.53 


0.6 


0.23 


1.2 


7.8 


0.35 


1.29 


1.27 


N2a 


505.3 


1.44 


6.44 


0.2 


0.15 


0.6 


8.2 


3.73 


0.73 


0.74 


N2b 


1.1 X 10'^ 


89.6 


14.5 


13.2 


0.84 


40.7 


62.5 


11.1 


2.80 


2.49 


N2c 


5.6 X 10^ 


8.75 


2.42 


1.4 


0.68 


0.7 


74.3 


19.4 


0.53 


0.54 


N2d 


2.9 X 103 


4.14 


6.48 


1.0 


1.43 


0.2 


74.7 


35.4 


0.34 


0.35 


N2e 


1.3 X 10^ 


9.00 


0.81 


1.2 


0.92 


1.5 


66.1 


18.2 


0.35 


0.36 


N3 


2.3 X 10^ 


3.79 


7.72 


1.7 


0.31 


2.0 


82.4 


11.8 


0.79 


0.79 


N4 


1.3 X 10^ 


3.29 


2.48 


1.4 


0.51 


0.8 


82.7 


11.8 


0.38 


0.39 


El 


111.6 


0.85 


0.77 


0.1 


0.11 


0.3 


2.2 


1.93 


0.39 


0.41 



to this question relies on magnetic reconnection. 

Magnetic reconnection occurs when two magnetic fluxes of opposite polarity encounter each other. In the presence 
of finite magnetic resistivity, the converging magnetic lines annihilate at the discontinuity surface and a current sheet 
forms there. In the standard Sweet-Parker (S-P) model the velocity at which the two converging fiuxes reconnect is 
given by Vrec ^ vaS^^'^^ where in this case gives the length of the reconnection layer (see Figure 15). Because S is 
large for Ohmic resistivity values (e.g., for the ISM, S ~ 10^^), the Sweet-Parker reconnection is very slow. In other 
words, since all the matter moving with the speed Vrec over the scale Lx must be ejected with the Alfven velocity 
through a thin slot (A), the disparity between the typical scales Lx and the outfiow thickness A, which in turn is 
determined by microphysics, i.e. the resistivity, makes the Sweet-Parker reconnection rate negligibly small. 

However, observations indicate that magnetic reconnection must be fast in some circumstances (e.g., solar fiares). 
Lazarian & Vishniac (1999) proposed a model for fast reconnection that is independent of the resistivity. The model 
appeals to the ubiquitous astrophysical turbulence as a universal trigger of fast reconnection. When turbulence is 
present within the current sheet, the outfiow region A gets determined by magnetic field wandering and therefore, 
becomes independent of the resistivity (see Figure 15, middle panel). It allows the formation of a thick volume filled 
with several reconnected small magnetic fluctuations which make the reconnection fast. This model was successfully 
tested numerically by Kowal et al. (2009) (see Figure 15, bottom panel). This challenges the well-rooted concept of 
magnetic field frozenness for the case of turbulent fiuios and provides an interesting way of removing magnetic fiux 
out of astrophysical fiows, e.g. star formation regions (Lazarian 2005; Santos-Lima et al. 2010; 2012; de Gouveia Dal 
Pino et al. 2011), accretion disks, or the solar dynamo. 

It must be remarked that numerical effects are always a concern when dealing with numerical simulations involving 
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Fig. 14. — The same as in Figure but showing the time evolution of the magnetic field lines of the turbulent model N2b from t = 1 
to 90 Myr. We clearly see how turbulence diffusively transports the magnetic field lines as time evolves. At t= 90 Myr the lines that were 
not transported to the outskirts of the cloud are advected to the center by the collapsing gas causing a local enhancement in the magnetic 
field intensity. 

reconnection and magnetic field diffusion. However, the high resolution numerical tests of magnetic reconnection 
performed by Kowal et al. (2009), showed that in the presence of turbulence the local non-linear enhancements of 
resistivity are not important. This is a confirmation that the turbulent reconnection diffusion that we observe in 
our simulations (see also Santos-Lima et al. 2010; 2012; de Gouveia Dal Pino et al. 2012) is a real effect and not a 
numerical artifact. Analytical studies summarized in Eyink et al. (2011) also support the notion that magnetic fields 
are generically not frozen-in when conductive fluids are turbulent. In view of these studies, we can conclude that the 
concept of reconnection diffusion looks very natural and ubiquitous. 
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